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O ■ Abstract 

A second order classical perturbation theory is developed and applied to elastic atom corrugated 
surface scattering. The resulting theory accounts for experimentally observed asymmetry in the 
C^ \ final angular distributions. These include qualitative features, such as reduction of the asymmetry 

I ' with increased incidence energy as well as asymmetry in the location of the rainbow peaks with 
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^ potentials. Analytic expressions for the angular distribution are derived for the exponential repul- 

sive and Morse potential models. The theory is implemented numerically to a simplified model of 
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'nI" ' the scattering of an Ar atom from a LiF(lOO) surface. 
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I. INTRODUCTION 



The scattering of atoms from surfaces has been measured extensively during the past fifty 
years jj-y]. The measured angular distribution for heavy atoms, when quantum diffraction 
effects may be neglected, is characterized by a few prominent, qualitative features. Rainbows 
appear at sub-specular and super-specular angles in the form of bell shaped maxima in the 
angular distribution p-LZj. Typica 
that of the super-specular peak 
distance between the rainbow peaks as the incident energy of the atom is increased 10|, lUl 
or as the angle of incidence (measured with respect to the vertical) increases |9|, |llj • Recent 



^ 



the intensity of the sub-specular peak is larger than 



10|. Other features include a reduction of the angular 



reviews of rainbow scattering from surfaces may be found in Refs. 12 



lal- 



An early model which gave a qualitative explanation for the asymmetry in the angular 
distribution was that of a hard wall corrugated potential 



14 



17l |. This model may be solved 



analytically. It provides a simple explanation for the rainbows - one readily finds that 
they originate from the inflection points of the corrugation. Perhaps more subtle but of 
not less interest is that the hard wall corrugated model also provides an explanation for 
the asymmetry in the angular distribution. One finds that the potential which faces the 
incoming particle leads to the sub-specular peak [l7[. Its intensity is higher just as the 
intensity of the rain that one feels is larger when one runs into the rain direction rather than 
away from it. 

The hard wall model also provides a partial explanation for the energy dependence of 
the distance between the rainbow angles. If one adds a shallow attractive square well which 
precedes the wall |lj, [l6|, |l7| , one finds that the distance between the rainbow peaks decreases 
as the incident energy increases. At low incident energies, one "feels" the shallow well and 
due to its refractive effect on the straight line trajectories, it increases the distance between 
the rainbow angles. As the energy is increased, the refraction decreases and one reaches the 
repulsive hard wall limit, in which the rainbow angles are energy independent. The energy 
dependence of the angular distance between the rainbows is thus a sensitive measure of the 
characteristics of the physisorption well. 

However, the hard wall class of models is deficient in a number of respects. For example, 
the distance between the rainbow angles is independent of the angle of incidence. The 
rainbow angles are symmetrically spaced about the specular angle and this is not always so 



[qI, |10|. Moreover, it is not very realistic, since the interaction between atoms and surfaces 
really is rather "soft" especially when dealing with rare gas incident atoms. 

With this in mind, we have developed in recent years a classical theory of atom surface 
scattering which is based on a perturbation theory in which the corrugation height is con- 



sidered to be the small parameter 
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2l|. In previous work, we developed this theory 



using a perturbation expansion which is valid to first order in the corrugation height. The 
first order theory correctly predicts the incident energy and incident angle dependence of the 
angular distribution. The distance between the rainbow angles decreases when the incoming 
atom has some time to traverse in the horizontal direction, smearing out the effect of the 
corrugation. Therefore, when the angle of incidence is large, the horizontal velocity is rela- 
tively fast and the distance between the rainbow angles is small. Similarly, at a fixed angle 
of incidence, increasing the energy implies also an increase of the horizontal momentum and 
this leads to a smaller angular distance between the rainbows. 

Similar to the hard wall model, if the potential is purely repulsive, then the rainbow angles 
are energy independent, the energy dependence arises only when the potential includes a 



physisorption well [13|. As in the hard wall model, the rainbows are symmetrically placed 
about the specular angle. The first order perturbation theory does not account for the 
asymmetry in the angular scattering unless one imposes an asymmetric corrugation potential 
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The topic of this paper is to show how a second order perturbation theory, applied to 
elastic atom-surface scattering, accounts for the asymmetry in the angular distribution. In 
principle, the second order theory calls for the solution of a second order in time equation of 
motion which is characterized by a time dependent harmonic frequency and external force. 
Such an equation is rather difficult to solve analytically except in special cases. However, if 
one uses the fact that energy is conserved during the scattering, one may replace the second 
order in time equation of motion with a first order in time equation which is readily solved 
analytically. Using this strategy we derive in Section II the second order perturbation theory 
expression for the final momenta of the particle. 

In Section III we use these results to derive the angular distribution expanded up to 
second order in the corrugation height. We find that the second order contribution leads 
to the correct asymmetry in the angular distribution. Sub-specular final angles have higher 
intensity than super-specular angles, due to the same qualitative effect already understood 



from the hard wall model. The perturbation theory also provides an incident angle and 
energy dependence for the asymmetry, which decreases as the energy is increased or as the 
angle of incidence increases with respect to the vertical direction. Moreover, the theory 
accounts for asymmetry in the location of the rainbow angles with respect to the specular 
scattering angle. 

In Section IV we apply the theory to an exponentially repulsive potential and a Morse 
potential model with a sinusoidal corrugation function. Analytic expressions are derived for 
the angular distributions in both cases. As in the first order theory, the purely repulsive po- 
tential gives an energy independent angular distribution. The energy dependence is directly 
related to the physisorption well, which is of course included in the Morse potential model. 
Some numerical examples are provided for a simplified model of the scattering of an Ar atom 
from a LiF(lOO) surface. We end with a discussion of the results, noting for example that 
one may also observe asymmetry in the opposite direction, that is that the super-specular 
scattering angles are more probable than the sub-specular ones [HI]. We also speculate that 
the second order perturbation expansion should be useful within a semiclassical description 
of the scattering. 



II. PERTURBATION THEORY FOR THE MOMENTA 

In this paper, we limit ourselves to a model of in-plane scattering from a frozen surface. 
Generalization to the full three dimensional dynamics, as well as inclusion of surface phonons 
is straightforward but leads to somewhat more complicated expressions [l3| . We thus assume 
that the scattering event takes place in the vertical (z) and horizontal (x) configuration space. 
A "standard" model used for the description of the scattering is based on the assumption 
that the potential of interaction depends on the instantaneous distance from the surface 
and so has the generic form V {x,z) = V {z — h {x)) where h (x) is the small periodic (with 
lattice length /) corrugation height of the surface. The potential vanishes when the particle 
is sufficiently distant from the surface. This potential is then expanded to first order in the 
corrugation height, that is: 

V{z,x) = V{z)-V'{z)h{x). (2.1) 

In principle, since we will develop a second order perturbation theory with respect to the 
corrugation, one should expand to include also the second order term, however, this does not 
create any fundamental differences in the results, only leads to a more complicated algebra, 
so that we will remain with the standard first order expansion for the potential. 

We will study the classical scattering for a particle with mass M and vertical and hori- 
zontal momenta Pz and Px respectively. The Hamiltonian governing the motion is thus: 

H = ^^^ + V{z,x). (2.2) 

The exact equations of motion governing the vertical and horizontal distances are 

Mzt + V'izt) - V"{zt)h{xt) = (2.3) 

Mxt-V'{zt)h'{xt) = 0. (2.4) 

where the primes denote derivatives with respect to the argument and the dots, time dif- 
ferentiation. The particle is assumed to be initiated at the time —to with initial vertical 
(negative) momentum p^- and (positive) horizontal momentum p^... The zero-th order mo- 
tion (expansion to order h^) is decoupled, the vertical motion is governed by the vertical 
Hamiltonian 

2 

H^o = ^ + Viz) (2.5) 



and the horizontal motion is that of a free particle 

HxQ = %. (2.6) 

with constant velocity Vx = Pxi/M. In the zero-th order motion, the particle impacts the 
surface at time t = and then leaves the interaction region by the time to which is taken to 
be sufficiently large to assure that the scattering event is over. Finally, we take the limit of 
to -^ oo. 

We then expand the horizontal and vertical motions to second order in the corrugation 
height: 

xt = xo^t + xi^t + X2,t + O [h (x) ) (2.7) 

Px,t = PxO,t + Pxl,t + Px2,t + O [h (x)^) (2.8) 

zt = zo^t + ^i,t + Z2,t + 0{h (xf) (2.9) 

Pz,t = PzO,t + Pzl,t + Pz2,t + [h (xf) . (2.10) 

The zero-th order vertical solution obeys the equation of motion: 

Mzo,t + V'{zo,t)=0 (2.11) 

and the first order correction to the vertical motion obeys the equation of motion: 

Mzi^t + V"izo,t)zi,t - V"izo,tMxo,t) = 0. (2.12) 

In the horizontal direction the motion is to zero-th order that of a free particle (parallel 
momentum conservation) such that: 

xo^t = Xo_to + j^ (^ + ^o) = ^O'O + j^t. (2.13) 

The Jacobian of the transformation between the initial value of the horizontal coordinate 
and its value upon impact Xo,o, is unity. The first order correction to the horizontal motion 
is determined by 

Mxi^t - V'izo,t)h'{xo,t) = (2.14) 

while the second order equation of motion for the horizontal coordinate is: 

Mx2,t - V'{zo^t)h" {xo,t) xi,t - V"{zo,t)zi,th' {xo,t) = 0. (2.15) 



One then readily finds that the horizontal momenta up to second order are expressed in 
terms of the lower order solutions as: 

P.i,t= [ dtV'izo,t)h'ixo,t) (2.16) 

J -to 

and 

P.2,t= [ dt[V'izo,t)h" ixo,t)xi,t + V"izo,t)z,,th' ixo,t)] . (2.17) 

J -to 

The corrugation function h (x) is periodic with period /. In the following, we will employ 
the simplest possible periodic corrugation function 

h{x) = hsm('^^ (2.18) 



but here too, we note that it is straightforward. 



harmonic expansion for the corrugation function 13|, ll9|]. Using the symmetry of the motion 



only increasingly complex to use a higher 



along the vertical direction {V'{zQt) is symmetric with respect to the time) we find that the 



first order contribution to the final horizontal momentum is 
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Pxij = Pz,K COS ( — -^ ) (2.19) 



with 

K = / dtV'{z^^t) cos {u^t) , (2.20) 

IPz, J-oo 

the horizontal frequency is defined as: 

--T^ (2.21) 

and we have taken the limit of to -^ oo. The dimensionless quantity K is termed the rainbow 
shift angle, as will also become evident from the expression for the angular distribution, as 
shown in the next Section. This first order contribution to the horizontal coordinate is then 
seen to take the form |l8| : 



Xi,t = COS 
where 



(^)/ dt'F^t')- sin l^^y^ dt'FAt') (2.22) 



27r r 
F,(t) = /i— / dt'V\zo,t')cos{coJ), (2.23) 

F^(t) = h^f dt'V {zo,t') sin {oo J). (2.24) 



The first order contribution to the vertical momentum necessitates the solution of Eq. 
12.121 This equation is equivalent in form to that of a forced oscillator with a time dependent 
frequency, which is rather difficult to solve analytically. We note however that to first order 
in the corrugation height, energy conservation implies that for any time 

= P-P-M +PmPzi,t ^ v\zoM,t - V'{zo,t)hixo,t). (2.25) 

This relation then determines the first order contribution to the vertical motion p^i^t in terms 
of the first order contribution to the horizontal motion: 

Pzij = -^Pxi,f- (2.26) 

Pzi 

Not less important is to note that the first order energy conservation relation 12.251 also 
leads to a first order in time equation of motion for the first order contribution to the vertical 
motion: 

zi,t = ^^zi,t + ^ r dth {xo,t) TX{zo,t) (2.27) 

PzO,t PzO,t J -to dt 

whose solution is readily written as: 



/■* 1 /■*' d - 

zi,t=P.o,t dt'^ dt"h{xo,t")-^V'{zo,t"). (2.28) 

J-ta PzQ.t' J -to "t 



dt'^— / 

-to PzO,t' J -to 

By taking the time derivative of Eq. 12.281 one may see that this is a specific solution to 
the second order in time equation of motion 12.121 for the first order correction. It is this 
observation that allows us to obtain closed form expressions for the final momenta, up to 
second order. Using the specific sinusoidal form for the corrugation (Eq. 12.181) Eq. 12.281 may 
be rewritten as 

z,,t = sin (^^) Pzo,t £ dt'G, (f) + cos f^:^\ p,o,t £ dt'Gs (f) (2.29) 



G^t) = -^ [ dt'^^^^^ cos (u J), (2.30) 



where we used the notation: 

PzO,t J-oo dt 

h_ /•* ,^,dV\zo,t' 

PzO,t 

With these preliminaries, using Eqs. I2.17[l?!^ and l?!^ we find after some manipulation 
that the second order contribution to the final horizontal momentum simplifies to: 

Px2,f = PxiKcc (2.32) 



Gs{t) = ^^ I dt' ^^ sin {io J'). (2.31) 



with 

K,^ = ^^ r dt [cos [u^t) ^^^^^ /* dt'Gs if) - ^V'izo,) sin (uj) f dt'F, [f 

(2.33) 

The second order contribution to the final vertical momentum is then found through energy 
conservation to be: 

P.2J = ^-V..j + %±^. (2.34) 

Pzi ^Pzi 

Eqs. 12. 32112. 341 are the central results of this Section. An explicit solution for the second 
order contribution to the final momenta and its dependence on the point of impact on the 
surface, has been derived. In principle, one can follow the same methodology to obtain 
all order contributions to the final momenta, however the complexity increases accordingly. 
One may also put in higher harmonics into the corrugation function, however, as we shall 
see in the next Section, the second order perturbation theory suffices to provide a qualitative 
explanation for the observed angular distributions, their asymmetry and rainbow structure, 
and angle of incidence and energy dependence. 



III. PERTURBATION THEORY FOR THE FINAL ANGULAR DISTRIBUTION 

The (negative) angle of incidence with respect to the vertical vector to the surface is by 
definition 



^. = tan-M^ . (3.1) 



The final angular distribution is (with x 



2TTX0fi - 



P,.,).^f«(«,-.a„-.(£^)) (3.2) 

where pxf {x) and pzf {x) are the final momenta as determined from Hamilton's equations 
of motion. Using the second order expansion results from the previous section, this may be 
rewritten as: 

with 

j/^N _ Vxi+PxlJ+Px2,f 



Pzi - PzlJ - Pz2,f 

K cos {x) + tan 6i{l + Kcc) 
1 — tan 9iK cos (x) — Kcc f an^ 9i 2 coT^ e^ 



(3.4) 



Using the notation y = K cos x we note that imposing that the argument of the Dirac 
"delta" function in Eq. 13.31 vanishes leads to a quadratic equation in y whose solutions are 
denoted as y*: 



, _ 1 ± y/l + 2 [tan {Oj + Of) - tan Oj] [tan {Oj + Of) + Kcc tan Oj] 
^^~ tan(% + ^i)-tan^i ^ ' ^ 

The condition that —K < y* < K and noting that the perturbation theory implies that the 
magnitudes of the rainbow shift angle should be small, typically such that \K\ < 1, implies 
that only the minus sign solution is physical. To gain some further insight, and staying 
consistent within the second order perturbation theory we may expand this solution: 

y* = yl + y; + 0{h') (3.6) 

where the subscript denotes first and second order in {9f + 6i) respectively, to find: 

yl = -tan(% + ^i), (3.7) 

^ / tan2(^. + ^i)\ , , 

yl = - tan Oi Kcc + H ■ (3.8) 
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To obtain the angular distribution it is necessary to determine the derivative of the 
deflection function at the points y* ,which may be written as: 

-— = —Ksmx-—. (3.9) 

ax ay 

Rainbows are found when the derivative of the deflection function vanishes that is when 

/' (x) = 0, provided that j- does not diverge at these points. One flnds that: 



dy 



df _ l + ^ + ytane,{l + K, 

dy 



cos^ 9i ll — tan 6iy — K^c tan^ 9i — j^ 
so that the rainbows occur either when 



2 



(3.10) 



sins = (3.11) 

or f = 0: 

dy 

l + yfltan^i(l + ifec) + Y = (3.12) 

where the subscript reminds us that this is the condition for the rainbow angles. For the 
condition sin (x) = we have two solutions, x = or a; = tt. In the second case, we have a 
quadratic equation with solutions 



yR,± = - tan^, (1 + K,,) ± y^tan^ 6, (1 + K,,)^ - 2. (3.13) 

These solutions will only be valid if 

tan^ 9i {1 + Kccf > 2 (3.14) 

and this implies collisions which are close to grazing collisions, or more speciflcally, \6i\ > 60°. 
The angular distribution is then 

cos^ e,[l- tan Oiy* - tan^ OiKcc - j^-rf 
P (^/) = ^ , , '^ • (3.15) 



vr cos^ 9 



l + rtan^,(l + i^ec) + ^ V^ 



*2 

y 



2 

To make further sense of this deceptively simple but in fact quite complicated expression 
we resort to perturbation theory. Expanding up to second order using the second order 
expansion of y* as in Eqs. 13.71 and 13.81 the angular distribution simplifles to 



, , 1 + tang, tan (Of + g,) + K,,ia.n^ Oi + "'"" \'^''> (tan^ g, - 1) , , 

P [Of) ~ ^-^ ^ ^ -. 3.16 

-K cos'' [9 f + ei)^K^-y*' 

11 



Rainbows occur whenever y* = ±K, but the strength of the divergence depends on the 
sign of y* due to the hnear term in the numerator. This gives rise to the asymmetry of the 
distribution. The rainbow condition imphes to second order that 

tan(% + ^i) = ±K -tanOi (Kcc+ — ) (3.17) 

showing that the location of the rainbows is no longer symmetric about the specular angle. 
Using the minus sign for ±1^1 corresponds to a rainbow angle which is smaller than the 
specular angle. Since the linear term in the intensity goes as tan 6i tan {6f + 6i) and 6i is 
negative, we get that tan^jtan(0j + 6i) > and the amplitude for the smaller scattering 
angle is increased relative to the amplitude for the larger scattering angle. If, as is typically 
the case, the rainbow shift angle i^ > 0, this implies (through the relation y = Kcos{x)) that 
X = 0. When considering the corrugation potential —V {z) hsin (x) this implies that when 
the particle hits the part of the potential which points towards the particle, the scattering 
intensity is higher. This is also qualitatively consistent with rainbow scattering from a hard 
wall corrugated potential. 

It is also worthwhile noting that when limiting the derivation to only first order pertur- 
bation theory (denoted by the subscript) then the angular distribution takes the symmetric 
form [lal,^: 

P (9f) = / (3.18) 

■^ TT cos2 {Of + 6,) y/K^ - tan^ {Of + Oi) 

demonstrating clearly that the asymmetry observed in the angular distribution is a second 
order effect. Finally it is instructive to compare the second order result Eq. 13.161 with 
the angular distribution found for scattering from a hard wall potential with the same sine 
corrugation function as in Eq. 12.181 



f 1 + tan f ^) tan {9 
Ph^ {Of) = ^ \ ^ ' = (3.19) 



7rcos2(^^j^irL-4tan2 

where 

Kn^ = ^. (3.20) 

The hard wall potential exhibits rainbow angles that are symmetrically placed around 
the specular angle. The asymmetry in the distribution which comes from the second term 
in the numerator is half as large as the asymmetric part in Eq. 13.161 Here one notes though 

12 



that the hard wall model is qualitatively different from the potential used in the present 
model. In the hard wall model, the location of the hard wall varies with the corrugation 
while the hard wall limit of the model we have been using (Eq. 12. ip leads to a hard wall 
whose location is independent of the horizontal coordinate. However, the qualitative feature 
of a larger scattering amplitude originating from the wall facing the incident particle is the 
same in both cases. 
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IV. ANALYTICAL MODELS AND NUMERICAL RESULTS 



A. Repulsive exponential potential 



Perhaps the simplest soft potential is obtained by replacing the hard wall with a purely 
repulsive exponential potential 



Ve{z) = Veexp{-az) 



(4.1) 



where the inverse length a is the stiffness parameter of the exponential potential. Vg has 
the dimensions of energy and expresses the " strength" of the exponential interaction. The 
trajectory for the vertical motion at energy Ez is known analytically 13|: 



exp (azf) 



A 
2K 



with 



n' 



[1 + cosh (nt)] 



2a^Ez 



M 



(4.2) 



(4.3) 



and Ez is the incident energy in the vertical direction. The rainbow angle shift (see Eq. 
I2.20p for this model is: 

(4.4) 



K. = K, 



TlQ 



hw , 



sinh (yrfi) 



with 



il = -— = — - tanb'j 



(4.5) 



and in the hard wall limit, that is when the stiffness parameter a ^ oo, K^ ^ K^y^. 

It is a matter of some algebra, based on the integral (for further manipulations see Section 
IVofRef. Iiol) 



dt- 



cos 



(Qt) 



and the identity: 



to find that: 



27rsinh($(]) 
cos $ + cosh (t)] sin $ sinh (vrO) 



[1 + cosh (t)] _ d / 2 cosht + 2 
sinh^ (t) dt \ sinh t 



K. 



Kl 



1 



■nVt cosh {VLti^ 



(4.6) 



(4.7) 



(4.8) 



tan^ 9i \ cos^ 9i sinh {fl-n) 
This implies, that for the repulsive exponential potential, also to second order in per- 
turbation theory, the angular distribution is independent of the incident energy. It only 



14 



depends on the angle of incidence. The angular distribution, to second order in the pertur- 
bation theory is given by Eq. 13.161 In the hard wall limit one has that the second order 
parameter K^ce -^ -Kl^. 

B. Morse potential model 

The dependence of the angular distribution on the incident energy comes from the ex- 
istence of a shallow physisorbed well in the potential of interaction of the atom with the 
surface. This is well modelled by the Morse potential 

Vm {z) = K) [(exp {-az) - if - l] (4.9) 

which has a physisorption well depth Vq. In this case, the rainbow shift parameter (Eq. 

[2:20]) is 

7rficosh(i7$) 

Km = ^/^^-^-r/^ (4.10) 

smh ynil) 

with 



and n as defined in Eq. 14.51 For a fixed angle of incidence, the angle $ decreases from vr to 
7r/2 as the energy is increased, causing the rainbow shift parameter to decrease accordingly. 
The attractive well leads to a narrowing of the distance between the rainbow angles, as the 
incident energy is increased. 

The trajectory for the Morse potential at an incident energy Ez is known analytically 

COS $ 

exp {azt) = 5— [cosh (fit) + cos $] , (4.12) 

sin $ 

with the frequency fl as given in Eq. 14.31 We then use the same integrals as for the 
exponential repulsive potential and the identity: 

1 _ a^ d /^fitsinh(fit) -cosh(fit) [l + cos2<I>] -2cos<I)\ 



pIq^, M^n^d{nt) V sinh(fit) 

to find that: 

K^^ = KIj ( 2IJT + .211 {^^ t^^h {^^) ^ ^^ ^°th (^^)) ~ t^^h {^^) ^ ^os $ sin $ J 

ytan |c/j| sin |c7j| J 

(4.14) 
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FIG. 1: (color online) The second order coefficient Kcc for the Morse potential model given by 
Eq. 14.141 is plotted as a function of the incident energy (in meV), for two incident angles, —45° 
and -30°. 

In the limit that the stiffness parameter a — )■ oo this reduces to the previous result, that 
is Kcc = ~Khw ■'■^ ^^^ limit of a purely repulsive exponential potential ($ — )■ 0) we regain 
the exponential model result. In the high energy limit, that is when $ — t- 7r/2 we have that 

1 n Q 



lim Kr 



K 



M 



tan \6i 



2 sin' % 



coth 



while in the low energy limit we have that $ — t- vr so that 

1 



lim Kr 



K 



M 



tan' \dj\ 



(4.15) 



(4.16) 



and this reduces to —K^,^ in the hard wall limit (a -^ oo). We thus find that the second 
order coefficient is positive at low energies and becomes negative as the energy is increased. 
Its magnitude reaches a maximum and then decreases with increasing energy. This leads to 
a reduction of the asymmetry in the angular distribution as the energy is further increased. 

C. Numerical examples 



The scattering of an Ar atom from the LiF(100 
both experimentally [10| as well as theoretically |l3| 



surface has been studied in some detail, 
. The experimental angular distribution. 
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FIG. 2: (color online) The second order coefficient K^c for the Morse potential given by Eq. 14.141 
is plotted as a function of the incident angle for two incident energies, 300 meV and 700 meV. 



measured with a fixed angle of 90 degrees between the incident beam and the detector showed 
a number of distinct qualitative features. The distance between the rainbow peaks decreased 
as a function of increasing incident energy. The asymmetry in the angular distribution was 
such that the intensity of the rainbow peak was higher for final angles which were less 
than 7r/4, however the asymmetry decreased as the incident energy decreased. All of these 
features are accounted for within the present second order perturbation theory. 

To demonstrate this we employ parameters used previously to fit the experimental results. 
In particular, the following values have been used for the Morse potential model 131] : h = 
0.25 a. n.,l = 4A,al = 3 and Vq = 88 meV. 

As has been previously stated, the second order coefficient, Kcc, plays a key role in the 
asymmetry of the angular distributions. The dependence of this coefficient on the incident 
energy and incident ang le is shown in FigO In the range of energies (315-705 meV) probed 
by the experiment |lO| the magnitude of the coefficient is a decreasing function of the 
incident energy. At very low energies, the asymmetry of the angular distribution is expected 
to be very important since this second order coefficient becomes relatively large. When the 
incident energy is very low, the approaching atom "feels" the corrugation for a longer period 
of time thus distinguishing between the case that the atom approaches the downhill or uphill 
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part of the corrugation. Given the large magnitude of the coefficient one should expect that 
in this low energy limit the perturbation theory will not be accurate. As the incident energy 
increases, K^c becomes small, approaching the purely repulsive model result, the asymmetry 
in the angular distribution is reduced and the perturbation theory result should be rather 
accurate. 

The variation of the second order coefficient with respect to the angle of incidence is 
plotted at two different incident energies in Figj2j These results have a number of interesting 
features. Firstly, at low angles of incidence and low energies the coefficient becomes positive. 
Secondly, the parabolic structure implies that similar values are obtained at different incident 
angles, indicating that the asymmetry is not necessarily a monotonic function of the angle 
of incidence. Fourthly, as may also be discerned from Eq. 14.141 when the angle of incidence 
tends to 7r/2 (grazing angle) the coefficient vanishes. 

Another property which emerges from the second order perturbation theory is that in 
contrast to the hard wall model, the location of the rainbows is no longer symmetrically 
distributed about the specular angle. A measure of this "rainbow asymmetry" is obtained 
by considering the difference between the angular distance of the rainbow angles from the 
specular angle. In the symmetric case, this difference of course vanishes. In Fig. [3l using Eq. 
I3.17[ we plot the rainbow asymmetry (in degrees) obtained by subtracting the distance of 
the subspecular rainbow peak from the specular angle from the distance of the superspecular 
rainbow angle from the specular angle. From this figure we note that, depending on the 
angle of incidence, the rainbow asymmetry can change sign. In addition the dependence of 
the asymmetry on the incidence energy is not necessarily monotonic. 

When considering the angular distribution, one distinguishes between two different ex- 
periments. In one class, the angle of incidence is kept fixed and the detector is moved to 
measure the final angular dependence of the outcoming fiux. In a different (easier) experi- 



mental setup, as used in the measurements of Ref. [10|, the angle between the incident and 
final beam is kept fixed and only the angle of incidence is varied. The results for the angular 
distribution for the former case are shown in Fig. H] for two different incidence energies. 
The singularities associated with the rainbow peaks have been smoothed by approximating 
the step function with a hyperbolic tangent function. As the hyperbolic tangent function 
tends to the step function the rainbow peaks become higher but it becomes more difficult to 
resolve their magnitude numerically. The present depiction suffices to show the asymmetry 
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FIG. 3: (color online) Rainbow asymmetry evaluated from Eq. 13. 171 for the Morse potential model 
and plotted as a function of the incident energy and three incident angles, -15°, -25° and -45° (see 
text). 



in the distribution. One notes that the asymmetry decreases significantly as the energy is 
increased from 300 meV to 700 meV. The fixed (7r/2) angular distribution, plotted as a 
function of the exit angle {Of = 7i/2 + 6i) is shown in Fig. [5l Qualitatively, the results are 
similar to those shown in Fig. HJ however the distance between the rainbow peaks becomes 
smaller. 



19 



P(e,) 





^1 


\ E. = 700meV J 


A 




E. = 300meV 



10 



20 



30 



40 



50 



60 



FIG. 4: (color online) The final angular distribution as given by Eq. 13.151 is plotted for two 
incident energies, 300 meV and 700 nieV (covering the experimental range) and for an incident 
angle of - 45°. 




FIG. 5: (color online) The final angular distributions as given by Eq. 13.151 but with a fixed angle 
of 7r/2 between incident and outgoing beams is plotted for two incident energies, 300 meV and 700 
meV, covering the experimental range. 
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V. DISCUSSION 

A second order perturbation theory with respect to the corrugation height has been 
developed for the elastic scattering of atoms from a periodic corrugated surface. The second 
order theory correctly accounts for experimentally observed asymmetry in the measured 
angular distributions. Expressions have been derived for the energy and incident angle 
dependence of the angular distributions and their asymmetry. In contrast to the hard wall 
models, the second order theory provides also the dependence of the asymmetry on both 
angle of incidence and energy of the particle. Analytical expressions for this dependence 
were derived for a purely repulsive exponential model potential as well as a Morse potential 
which exhibits the characteristic physisorption well felt by the incoming atom. Numerical 
results were shown for parameter values which fit qualitatively the scattering of Ar from a 
LiF(lOO) surface. In contrast to the corrugated hard wall model, the second order theory 
also accounts for asymmetry in the location of the rainbow angles. 

In this context it should be noted though that the asymmetry is not always such that 
the subspecular rainbow peak is the preferred one. The opposite is found for the scattering 



of Ar from an H covered Tungsten surface ll| . This indicates for example that the elastic 



theory presented in this paper, may not always be sufficient. For example, phonon friction 
which is larger in the vertical direction than in the horizontal direction will tend to shift the 
final angular distribution towards superspecular peaks. The three dimensional structure of 



the surface can also affect the asymmetry in the distribution [13|. 

As already noted, the present theory was limited to elastic scattering. It can be further 
expanded to include the interaction of the particle with surface phonons. One should expect 
that a first order theory with respect to the coupling to the surface phonons should suffice, at 
least when considering the angular distribution. It is though possible but rather cumbersome 
to follow the methodology presented in this paper to also treat the coupling to the phonons 
to second order in the coupling constant. 



The classical first order perturbation theory was used in previous work [13|, |22| to study 
the sticking of atoms scattered from surfaces. Here too, one could employ the present 
second order theory to study sticking. It would be of interest to understand how much the 
asymmetry will change the sticking probabilities. 

The first order perturbation theory fails especially for grazing angles, where the change in 
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the horizontal moinentuin can no longer be considered as small with respect to the magnitude 
of the incident horizontal momentum. The second order perturbation theory should improve 
the theory but the extent is not clear. Detailed comparison with numerically exact classical 
mechanics simulations of the scattering would be helpful in this respect. 

Finally, we note that the first order perturbation theory has been used extensively within a 
semiclassical context 23|,|2J]. It should be of interest to see whether the present second order 
perturbation theory can be employed semiclassically, so that also the resulting semiclassical 
diffraction patterns will exhibit the correct asymmetry. 
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